A Search for Time Variation of the Fine Structure Constant 
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A method offering an order of magnitude sensitivity gain is described for using quasar spectra to 
investigate possible time or space variation in the fine structure constant a. Applying the technique 
to a sample of 30 absorption systems, spanning redshifts 0.5 < z < 1.6, obtained with the Keck I 
telescope, we derive limits on variations in a over a wide range of epochs. For the whole sample 
Aa/a = —1.1 ±0.4 x 10~ 5 . This deviation is dominated by measurements at z > 1, where Aa/a = 
— 1.9 ± 0.5 x 10~ 5 . For z < 1, Aa/a = —0.2 ± 0.4 x 10" 5 , consistent with other known constraints. 
Whilst these results are consistent with a time- varying a, further work is required to explore possible 
systematic errors in the data, although careful searches have so far not revealed any. 
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There are several theoretical motivations to search for 
space-time variations in the fine structure 'constant', a. 
Theories which attempt to unify gravity and other fun- 
damental forces may require the existence of additional 
compact space dimensions. Any cosmological evolution 
in the mean scale factor of these additional dimensions 
will manifest itself as a time-variation of our bare three- 
dimensional coupling 'constants', Alternatively, the- 
ories have been considered which introduce new scalar 
fields whose couplings with the Maxwell scalar F a i,F ab 
allow time-varying a 12] . The measurement of any vari- 
ation in a would clearly have profound implications for 
our understanding of fundamental physics. 

Spectroscopic observations of gas clouds seen in ab- 
sorption against background quasars can be used to 
search for time variation of a. Analyses involving op- 
tical spectroscopy of quasar absorbers have concentrated 
on the relativistic fine-structure splitting of alkali-type 
doublets; the separation between lines in one multiplet is 
proportional to a 2 , so small variations in the separation 
arc directly proportional to a. to a good approximation. 

Whilst the simplicity of that method is appealing, the 
relativistic effect causing the fine splitting is small, re- 
stricting the potential accuracy. We demonstrate below 
how a substantial sensitivity gain is achieved by compar- 
ing the wavelengths of lines from different species, and 
develop a new procedure, simultaneously analysing the 
Mgll 2796/2803 doublet and up to five Fell transitions, 
(Fell 2344, 2374, 2383, 2587, 2600A), from three differ- 
ent multiplets. These particular transitions are chosen 
because: (i) they are commonly seen in quasar absorp- 
tion systems, (ii) they fall into and span a suitable rest- 
wavelength range, (iii) an excellent database was avail- 
able ||, (iv) extremely precise laboratory wavelengths 
have been measured, (v) the large Fe and Mg nuclear 
charge difference yields a considerable sensitivity gain. 

We describe the details of the theoretical develop- 



ments in a separate paper JB| , here summarizing the main 
points. The energy equation for a transition from the 
ground state within a particular multiplet, observed at 
some redshift z, is given by 



E z = E a 
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where Z is the nuclear charge, L and S are the electron 
total orbital angular momentum and total spin respec- 
tively, and E c is the energy of the configuration centre. 
The term in the coefficient Q\ describes a relativistic cor- 
rection to E c for a given change in a, ao is the zero 
redshift value, and a z is the value at some redshift z. 
Re-arranged, this gives 
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Eq.(^|) is an extremely convenient formulation, the sec- 
ond and third terms contributing only if a deviates from 
the laboratory value. Accurate values for the relativis- 
tic coefficients, Qi,Ki and Ki, have been computed us- 
ing relativistic many-body calculations and experimental 
data. The coefficients and laboratory rest wavelengths 
are given in eq.(|J). For Fell, the relativistic coefficients, 
(Qi), are at least one order of magnitude larger than the 
spin-orbit coefficients, (Ki). The variation of the Fell 
transition frequencies with a is thus completely domi- 
nated by the Qi term. In Mgll, the relativistic correc- 
tions are small due to the smaller nuclear charge Z (see 
eq.(|^)), so whilst a change in a induces a relatively large 
change in the observed wavelengths of the Fell transi- 
tions, the change is small for Mgll. The relative shifts 
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are substantially greater than those for a single alkali 
doublet alone (such as Mgll), so Mgll acts as an 'an- 
chor' against which the larger Fell shifts are measured. 
A comparison of the observed wavelengths of light and 
heavy atoms thus provides a dramatic increase in sensi- 
tivity compared to analyses of alkali doublets alone. 

Since it is already clear from previous observational 
constraints that any change in a will be very small Q, 
0, ||, it is vital that E z=0 is known accurately enough. 
Indeed, the change in the frequency interval between 
Mgll 2796 and Fell 2383 induced by a fractional change 
Aa/a = 10~ 5 is, using eq.(||), 0.03 cm -1 . Thus, indepen- 
dent of the quality of the observations, the limiting accu- 
racy in a determination of Aa/a is ~ 10~ 5 for an uncer- 
tainty in the laboratory frequency of ~ 0.03 cm" 1 . This 
highlights the advantage of comparing light and heavy 
atoms. Previous analyses of alkali doublets, ||], @, have 
used frequencies of about the accuracy above, but have 
been restricted to placing limits of Aa/a w 10 -4 . 

Very precise laboratory spectra of the Mgll 2796 and 
Mgll 2803 lines have recently been obtained jl3| , in excel- 
lent agreement with previous accurate measurements of 
Mgll 2796 alone @, p[. Similarly precise Fell hollow- 
cathode spectra exist JL2[. Inserting these laboratory 
wavelengths and our Q and K coefficients into cq.(|^), 
we obtain the dependence of frequency on a for Mgll 
(top two equations in eq.(||) below) and Fell: 
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The astronomical data used for this analysis was ob- 
tained using the HIRES echelle spectrograph [H on 
the Keck I 10m telescope during three observing runs 
in 1994-1996. High-quality spectra of 25 quasars were 
obtained, in which intervening absorption systems at 
low/intermediate redshift have been identified containing 
Fell, Mgll, and other species. Full observational details 
are given in ref . . 

We now determine the relative positions of the Fell and 
Mgll lines and estimate Aa/a for each absorption sys- 
tem in the sample. Measuring each line (ie. Mgll 2796, 
Mgll 2803 and 5 Fell lines) independently is not op- 
timal, because the number of fitting parameters is not 
minimised, as discussed below. The procedure used is it- 
erative, where all available lines are fitted simultaneously 
with Voigt profiles, using VPFIT, a non- linear least- 
squares programme designed specifically for analysing 
quasar absorption spectra pE| . We minimise the total 
number of free parameters by linking physically related 



parameters: (i) the redshifts of the corresponding Fell 
and Mgll components are tied; (ii) the column densities 
N(FeII), N(MgII) for individual components can vary in- 
dependently of each other, but the velocity dispersion 
parameters 6(FeII), 6(MgII), (6 = \/2a), in correspond- 
ing components, are constrained by yj (24/56) b (Mgll) < 
fc(FeII)< fc(Mgll). We carry out 150 separate fits to 
each absorption system, varying a slightly each time, us- 
ing eq.(||) to compute the input rest-frame wavelengths. 
That procedure is performed 3 times; twice where the 
b parameters are related according to the two extremes 
above (ie. thermal or turbulent line-broadening) and a 
third time where all b parameters vary independently. 
The fitting procedure returns a value of Xmin which was 
computed as a function of Aa/a. We used the standard 
statistical procedure of estimating ler errors on Aa/a 
from y 2 ■ i 1. 

Several consistency checks are imposed before accept- 
ing a result. First, the Xmin f° r eac h individual spectral 
region being fitted must be statistically acceptable (ie. 
its reduced y 2 is ~ 1). It follows that Xmin f° r the fit 
as a whole is statistically consistent with the number of 
degrees of freedom for that fit. Second, we require sta- 
tistical consistency between the three separate analyses 
(for the three different b constraints). If any of the three 
Aa/a differ by more than la from either of the other 
two, the system is rejected. These criteria lead to the 
rejection of only about 1/10 of the sample, reflecting the 
good general statistical robustness of the procedure. The 
final Aa/a accepted corresponded to the lowest of the 3 
values of x 2 min - 

The best-fit values of Aa/a, and the la error bars, 
are plotted against redshift in Fig. 1. For z < 1 there 
is no departure from the present-day value. The scatter 
in the data is consistent with the individual error bars, 
ie. there is no evidence for any space or time-variations 
in a. However, at z > 1 the situation is less clear, with 
14 points giving x 2 = 34.9, 3 falling above zero and the 
rest below. Fig. 1 shows there is a dip in a relatively 
narrow redshift interval 0.9 < z < 1.2, where 9 out of the 
10 points lie below zero, and it is this dip which cause 
the overall z > 1 % 2 to be high. Assuming that this is a 
statistical fluctuation and that we have somehow under- 
estimated the errors, we can increase the error bars on 
each point by a constant amount, S, where S > such 
that the reduced x 2 = 1, *e. the error on the i th point 
becomes a\ = ai + S (where <7j is the value illustrated in 
Fig. 1). The results of this procedure, for the sample as 
a whole, and in two redshift ranges are: 

= -1.09 ±0.36 x 10~ 5 (0.6 < 2 < 1.6) 
= -0.17 ±0.39 x 10~ 5 (0.6 < 2 < 1.0) (4) 
^ = -1.88 ±0.53 x 10- 5 (1.0 < 2 < 1.6) 

The values of S are, in order, S = 0.46, 0.06, 0.53. The 
whole sample departs from zero at the 3.0<r level towards 
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smaller values of a. However, the z < 1 points alone show 
no significant trend, suggesting that there are no sig- 
nificant errors in the adopted laboratory Fell and Mgll 
wavelengths, or in the general procedure. However, at 
z > 1, there is a 3.5tr deviation. We have also exper- 
imented with other statistical methods which also give 
more conservative error estimates than a simple weighted 
mean (eg. a Bayesian method), and obtain results con- 
sistent with the above (to be reported elsewhere). 

What systematic (non-physical) effects could mimic- 
such an effect (either a general trend towards negative 
Aa/a for z > 1, or, the curious 'cusp' in the range 
0.9 < z < 1 .2)? The quality of the spectroscopic data 
does not deteriorate with redshift (this would be revealed 
by larger error bars at higher redshift) but in fact im- 
proves The unmodified error bars reflect uncertain- 
ties due to signal-to-noise, spectral resolution, velocity 
structure, the number of absorption lines fitted, and as- 
sume we have properly de-convolved the absorption fea- 
tures into their correct number of individual components. 
However this assumption may be wrong. There may be 
weak blended unresolved lines, even with data of the high 
quality we have here. For a specific absorption system, 
an interloper could indeed mimic a shift, eg. if it were 
blended with one of the Mgll 2796 or 2803 'anchors'. 
Nevertheless it is unlikely that this occurs preferentially 
for z > 1 and not for z < 1. Nor should random blend- 
ing create a negative Aa/a trend; it should merely in- 
crease the scatter about a mean of zero. An extensive 
search for a systematic blend (ie. a weak transition from 
some other species at the same redshift as the absorbing 
cloud) falling close to one of the Mg or Fe absorption 
lines, proved negative (details to be reported elsewhere) . 

Uncertainties in the Fe or Mg laboratory wavelengths 
cannot be responsible, since any errors would be an or- 
der of magnitude below the observed effect, and in any 
case, the z < 1 points reveal no offset. We also checked 
to see whether uncertainty in the instrumental resolution 
could introduce errors, and found the results insensitive 
to this. It is conceivable that certain absorption redshifts 
could be less reliable, due to the positioning of the Fell 
and Mgll lines with respect to the echelle order edges, 
where wavelength calibration may be worse. At higher 
wavelengths, the HIRES optical format is such that there 
are gaps in the wavelength coverage. Also, the ThAr cal- 
ibration spectrum has fewer lines per unit wavelength 
at higher wavelengths. These effects could mean that 
wavelength calibration is less reliable for higher redshift 
points than the lower redshift ones. Alternatively, per- 
haps there is a subtle bug in the wavelength calibration 
software. The z > I effect we have found is very small: 
a gradual drift of ~ 2.5 times the wavelength calibration 
residuals, over the range corresponding to the observed 
Fell and Mgll lines, could produce an effect of the signif- 
icance we have found. Previous analyses of similar qual- 
ity data have not demanded such high precision, so such 



errors could have gone unnoticed in other spectroscopic 
analyses. 

Can any of these effects produce the apparent trend 
we find? We have carried out extensive numerical ex- 
periments aimed specifically at testing for these poten- 
tial systematic problems. These involved carrying out an 
identical analysis on the wavelength calibration spectra 
(ThAr) as for the QSO spectra, measuring the emission 
line wavelengths and inferring a 'change' in a from shifts 
betwen the measured and laboratory wavelengths (de- 
tails to be reported elsewhere). The results were unam- 
biguous, in that any errors in the wavelength calibration 
across all spectral orders are so small that they could not 
contribute significantly to the apparent offset at z > 1, 
unless the literature values for the calibration spectral 
lines (ThAr) are substantially in error, with a gradual 
non-linear shift in the ThAr wavelengths. 

We consider all of the above potential errors to be 
rather unlikely. Since the effect is dominated by the 
0.9 < z < 1.2 points, and because a genuine physical 
effect confined to one specific epoch in the history of the 
universe does not seem at present to be well motivated 
by theoretical expectations, we presume the explanation 
is that there are additional undiscovered velocity compo- 
nents in the absorbing gas for those particular absorption 
sytems, even though all points have been subjected to the 
same analysis. Further observations of different species 
in these clouds could answer this. Our results should 
thus be regarded as stringent upper limits on any pos- 
sible time- variation rather than a positive detection of a 
change. 

We may compare our results with other recent values. 
Observations at z ~ 3 have yielded an upper limit 
|Aa/a| < 3.5 x 10~ 4 . Our analysis is at lower redshift, 
so a comparison requires choosing some (arbitrary) func- 
tional form of the evolution. At lower redshifts, a recent 
analysis of radio wavelength spectra of atomic hydrogen 
and molecular gas M gives limits |Aa/a| < 3 x 10~ 6 at 
two redshifts, z — 0.25,0.68. An additional recent con- 
straint of 3.5 ±5.5 x 10~ 6 comes from a comparison of 
hyperfine and optical redshifts. These limits are consis- 
tent with the results given in eq.@. The strongest ter- 
restrial constraint on the time evolution of a comes from 
the Oklo natural nuclear reactor Jl6). The Oklo event 
is estimated to have taken place around 1.8 billion years 
ago (corresponding to z w 0.1). We adopt a cosmological 
model with go = 1/2, A = and take the age of the uni- 
verse to be 13 x 10 9 years. The Oklo data has recently 
been re-examined [ fl7| , and upper limits on a change in 
a are -0.9 x 10~ 7 < (a okl ° ~ a )/a„ < 1.2 x 10~ 7 , as- 
suming the weak and strong couplings are unchanged. 
These bounds are clearly consistent with our results for 
z < 1. For the cosmological parameters quoted, if we 
adopt z = 0.8, the mid-point of our lower redshift range, 
cq.(^) implies 
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= -2.2 ±5.1 x 10" 
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yr 1 (0.6 < z < 1.0) (5) 

What other physical phenomena, other than time vari- 
ability of a, could give rise to the observational effect we 
report? The spacing of the Mgll and Fell isotopes is such 
that a significant change in the isotopic ratios could ex- 
plain the observations. However, the change would need 
to be substantial; for example, this would require most of 
the Mg in the universe at z > 1 to be in 26 Mg (the present 
epoch abundance is ~ 10%), and a physical mechanism 
found to convert almost all the 26 Mg into 24 Mg by the 
present epoch. If large scale magnetic fields exist, and 
the quasar light is polarised, these could potentially give 
rise to correlated apparent shifts in absorbers in neigh- 
bouring regions of the universe. However, for magnetic 
fields to be responsible for the global effect in a for z > 1 
(but not for z < 1), a sharp variation at z ~ 1 or some 
form of oscillatory variation would be required, both of 
which are hard to motivate. 

Some authors have suggested more exotic forms of evo- 
lution of the constants, including oscillations [|J. These 
could arise from new light bosons with mass m, pro- 
ducing periodic variations in the frequency of the radi- 
ation emitted at high z, with a modulation frequency 
~ vnT 1 ||. The creation of time-varying a by approxi- 
mate global symmetry ( j|]) allows oscillatory variations 
introduced by the decaying mean oscillations of a scalar 
field coupling to F a f,F ab . We note that as we move 
above z — 1.25, for a critical-density non- inflationary 
universe, we can encounter causally disconnected regions 
of the universe. The observational sensitivity achieved 
by these observations exceeds that of current microwave 
background observations and a larger dataset than ours 
may contain important new cosmological information. 

The work we have presented here demonstrates the 
possibilites for extending this type of study to incorpo- 
rate different species, other than Mg and Fe, spanning a 
wider redshift baseline, so equally impressive constraints 
should be obtainable at higher redshifts. Future analy- 
ses of other species will be hampered by the lack of ac- 
curate laboratory wavelengths. We hope that this paper 
will provide an impetus for new high-precision laboratory 
measurements. 
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FIG. 1. Plot of Aa/a vs z for 30 Fell/Mgll absorption 
systems. All the transitions given in equation 3 were used, 
where available. Plotted error bars were determined from 
Xmin ± 1 (but statistical results estimated using larger errors, 
as discussed in the text). 
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